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Abstract 

We investigate the relation between the phase space structure of Hamiltonian and non- 
Hamiltonian deterministic thermostats. We show that phase space structures governing reac- 
tion dynamics in Hamiltonian systems map to the same type of phase space structures for the 
non-Hamiltonian isokinetic equations of motion for the thermostatted Hamihonian. Our results 
estabhsh a framework for analyzing thermostat dynamics using concepts and methods developed 
in reaction rate theory. 
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I. INTRODUCTION 



Deterministic thermostats are widely used to simulate equilibrium physical systems de- 
scribed by ensembles other than microcanonical (constant energy and volume, (EA^)), such 



)lume, (l 



De- 



as constant temperature- volume (T, V) or temperature-pressure (T, p) 
terministic thermostats are typically obtained by augmenting the phase space variables of 
the physical system of interest with a set of additional variables whose role is to alter the 
standard Hamiltonian system dynamics in such a way that a suitable invariant measure in 
the system phase space is preserved. In the familiar Nose-Hoover (NH) thermostat , for 
example, the exact dynamics preserves both an extended energy Ti and a suitable invariant 
measure, ensuring that, provided the extended system dynamics is efectively ergodic on the 
timescale of the simulation, the physical system will sample its phase space according to the 
canonical (constant T) measure. 



non- 



Extended system thermostat dynamics can 
-Hamiltonian 114,1151, 
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Ul or 



20| . An important motivation for the 



brmulation of Hamiltonian deterministic thermostats such as the Nose-Poincare system 



111 is the possibility of using symplectic integration algorithms to compute trajectories 



21 



221. 



In this approach, an extended Hamiltonian is defined for the physical system plus thermo- 
stat variables which incorporates a coordinate-dependent time scaling of Poincare-Sundman 



type 



23 



24( 1 . Restricting the dynamics to a fixed value (zero) of the extended Hamiltonian 



results in the system variables sampling their phase space according to, for example, the 



canonical density 



111 (subject to the assumption of ergodicity.) The Hamiltonian version 



of the isokinetic thermostat is described in Section [TTl below. 

A fundamental question concerning deterministic thermostats has to do with the effective 
ergodicity of the dynamics on the timescale of the simulation. If the dynamics is not effec- 



tive 



|25|, 



y ergodic, then trajectory simulations will not generate the correct invariant measure 



26| . It has long been recognized, for example, that the dynamical system consisting of 
a single harmonic oscillator degree of freedom coupled to the NH thermostat variable is not 
ergodic l?! . A large amount of effort has been expended in atternpts to design thermostats 



exhibiting dynamics more ergodic than the basic NH system 



m 



29|. 



The question of ergodicity in thermostats is conceptually closely related to the problem 



2 



of statistical 

3, 



rates 



31 



versus nonstatistical behavior in the (classical) theory of unimolecular reaction 



32| . Broadly speaking, in this case one would like to know whether a molecule 



will behave according to a statistical model such as RRKM theory, or whether it will exhibit 



33, 



3|. Such 



significant deviations from such a theory, ascribable to nonstatistical dynamics 
'nonstatisticality', which can arise from a number of dynamical effects, is analogous to the 
failure of ergodicity in deterministic thermostats. 

In recent years there have been a number of theoretical and computational advances 
in the application of dynamical systems theory 35|, |36|, |37t| to study reaction dynamics 



and phase space structure in multimode models of molecular systems, and to probe the 
dynamical origins of nonstatistical behavior 38|, |39|] . The fundamental chemical concept of 
the transition state, defined as a surface of no return in phase space, has been successfully 



and rigorously generalized from t 
with > 3 degrees of freedom 



le well-established 2 degree of freedom case [40!] to systems 



39] . Moreover, dynamical indicators exist (determination 



of reactive phase space volume, behavior of the reactive fiux) to diagnose nonstatistical 
behavior. 

Despite their obvious potential relevance for the questions at issue, there has been rela- 
tively little work applying the powerful techniques from modern dynamical systems theory, 
in particular the theory of multidimensional Hamiltonian systems 

4^ 



36 



37 1 , to study the phase 



44|. There appears 



46, 



43] to the 



space structure of deterministic thermostats [2], [s], 0, 29. |41. 
to be considerable scope for application of these and other approaches [4| 
dynamics of deterministic thermostats. 

In the present paper we begin the development of a novel theoretical framework for the 
study of thermostat dynamics. Specifically, we describe how recently developed methods 
for the analysis of multimode Hamiltonian systems can be applied to investigate the phase 
space structure of the isokinetic thermostat [2|]. 

Although not as widely used as the Nose-Hoover thermostat and its many variants, the 
non-Hamiltonian version of the isokinetic thermostat has been developed and applied to 



several problems of chemical interest by Minary et al. 48 



49[. In this thermostat, the 



particle momenta are subject to a nonholonomic constraint that keeps the kinetic energy, 
hence temperature, constant. The resulting dynamics generates a canonical distribution 
in configuration space j^. A Hamiltonian version of the isokinetic thermostat was given 



by Dettmann 



and this Hamiltonian formulation (see also 



50 



51| ) is the point of 



3 



departure for our investigation. 

The Hamiltonian formulation of the isokinetic thermostat is presented in Section [III The 
non-Hamiltonian equations of motion for a Hamiltonian system subject to the isokinetic con- 
straint are shown to correspond to Hamiltonian dynamics at zero energy under an extended 
Hamiltonian whose potential is obtained from the physical potential by exponentiation. 
The extended Hamiltonian dynamics are therefore nonseparable and potentially chaotic (er- 
godic), even though the physical Hamiltonian might be separable. For the Hamiltonians we 
consider the physical potential exhibits a saddle of index one, as for the case of a bistable 
reaction profile coupled to one or more transverse confining modes. The bistable mode can 
play two distinct roles in the theory: it can either be interpreted as a reaction coordinate of 



48|. 



physical interest, or as a thermalizing thermostat mode 

Essential concepts concerning the phase space structure of multimode Hamiltonian sys- 
tems, especially the significance of Normally Hyperbolic Invariant Manifolds (NHIMs) and 
their role in the phase space structure and reaction dynamics of multimode molecular sys- 
tems with index one saddle, are briefly reviewed in Section IIIII 

In Section[IV]we show that the extended Hamiltonian dynamical system satisfies the same 
conditions satisfied by the physical Hamiltonian that give rise to the phase space structures 
discussed in Section IIIII We then show that these phase space structures exist for the non- 
Hamiltonian isokinetic equations of motion for the thermostatted physical Hamiltonian by 
an explicit mapping. Section |V] concludes. 

II. THE PHYSICAL HAMILTONIAN, THE EXTENDED HAMILTONIAN, AND 
NON-HAMILTONIAN ISOKINETIC THERMOSTAT 

We begin with a physical Hamiltonian of the standard form 

H{q,p) = ^-p' + ^q), (1) 



with Hamilton's equations given by: 



<i = ^=P- (2a) 
P=-^ = -^(^), (2b) 
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where {q,p) G M"" x are the physical coordinates and is the potential energy. Follow- 
ing Dettmann and Morriss [2|,]8l, we construct a Hamiltonian system having the property 
that trajectories on a fixed energy surface of the new Hamiltonian correspond to trajecto- 
ries of the physical Hamiltonian ([T]) which satisfy an isokinetic constraint in the physical 
coordinates. An extended Hamiltonian K, is defined as follows: 

/C(g,7r) = e-^*?i>B (3) 

where is 

7i<B = ie('»+i)*7r2-ie(^-i)*, (4) 

*B is an arbitrary parameter, and the relation between momentum variables p and tt is 
specified below. The value chosen for the parameter 03 defines a particular time scaling via 
factorization of /C; setting 03 = —1, for example, ensures that H<b has g-independent kinetic 
energy. (For simplicity we measure energies in units of ksT, thus keeping the value of T 
implicit.) 

The Hamiltonian ([3]) includes a time scaling factor e"*^*, and Hamilton's equations of 
motion for /C in physical time t are: 

die ^ 

q = +— = e TT (5a) 



die 



1 $ 2 , 1 -$ 

-en + -e 

2 2 



(5b) 



and are manifestly Q3-independent. 

To show that trajectories of the Hamiltonian system with /C = correspond to 
trajectories of the physical system ([T]) satisfying the isokinetic constraint, first note that the 
time derivative of T^sg along trajectories of ([5]) is given by 

• . dH^ . dH^B /„ X 

= Q -t: h TT — — (6a) 

oq on 

= "B^ge^nH^. (6b) 

This implies that trajectories of ([5]) satisfying Tisg = at t = 0, satisfy 7i«8 = for all t (for 
arbitrary ^). Using ([3]), this implies that these trajectories are also confined to the surface 
/C = for all time. 

The relationship between the Hamiltonian dynamics of ([5]) on /C = and isokinetic 
trajectories of ([1]) is made apparent by making the noncanonical transformation of variables 



q, (7a) 
7, H-> e-^^^V- (7b) 

This coordinate transformation is clearly invertible and is, in fact, a diffeomorphism (as 
differentiable as $). Applying ([7]) to Tifs gives: 

V-^^*(/-l) (8) 

from which we can immediately conclude that trajectories of the Hamiltonian system ([5]) on 
K. = ?i«8 = automatically satisfy the isokinetic condition 

= 1, (9) 

in the physical coordinates {q,p)- Substituting relation ([7]) into ([5]) we obtain equations of 
motion for {q,p): 

q = p (10a) 

p = ^(p2 + 1) + P(*g • q) = - ap (10b) 

where a = —^q-p and we have used the constraint p'^ = 1. Equations ( fTOj) are the isokinetic 
equations of motion for the thermostatted physical Hamiltonian ([1]) in physical time t, 
obtained via Gauss' principle of least constraint 

By design, the isokinetic dynamics ffTOj) generates a canonical distribution in the coordi- 
nates 0, ^- Minary et al. 48| have shown that the addition of thermalizing degrees of 
freedom to the physical Hamiltonian ([T]) can facilitate the attainment of the correct canon- 
ical distribution in g-space. If H describes a collection of uncoupled oscillators, addition 
of a bistable thermalizing degree of freedom renders the Hamiltonian dynamics under K, 
isomorphic to that of a reactive degree of freedom coupled to several bath modes, so that we 
can obtain useful insights into the thermostat dynamics using methods recently developed 
for multidimensional Hamiltonian systems. Alternatively, if H describes a reactive mode 
coupled to bath modes, then the /C dynamics is already in an appropriate form for the phase 
space analysis described in the next Section. 
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III. PHASE SPACE STRUCTURES ON A FIXED ENERGY SURFACE 



Our analysis of thermostat dynamics will be carried out in phase space, using the tools 



and framework for reaction type dynamics of Hamiltonian systems developed in 



38 



53 



54 



55 



57 



3, 



5d, 



58| . We will show in Section [IV] that these results apply both to the 
physical Hamiltonian system ([2]) and the extended Hamiltonian system ([5]). Here we give a 
brief summary of the setting and relevant results from these references. 

The starting point for identifying a region of phase space relevant to reaction is to locate 
an equilibrium point of Hamilton's equations, denoted {q*,p*), that is of saddle-centre-. . .- 
centre stability type. By this we mean that the matrix associated with the linearization of 
Hamilton's equations about this equilibrium point has two real eigenvalues of equal mag- 
nitude, with one positive and one negative, and n — 1 purely imaginary complex conjugate 
pairs of eigenvalues. We will assume that the purely imaginary eigenvalues satisfy a generic 
nonresonance condition in the sense that they are independent over the rational numbers 
(this is discussed in more detail in Section HVl) . 

We will assume that such an equilibrium point is present in the physical system ([2]) and 
we will show that the same type of equilibrium point exists for the extended Hamiltonian 
system ([5]) in Section [IVl However, the discussion in this section applies to any type of 
Hamiltonian system near the same type of equilibrium point. Without loss of generality we 
can assume that {q*,p*) is located at the origin, and we denote its energy by H{q*,p*) = h*. 

We will be concerned with geometrical structures in a neighborhood of phase space con- 
taining the saddle-centre-. . .-centre type equilibrium point. We emphasize this fact by de- 
noting the neighborhood by C; this region is to be chosen so that a new set of coordinates can 
be constructed (the normal form coordinates) in which the Hamiltonian can be expressed 
(the normal form Hamiltonian) such that it provides an integrable nonlinear approximation 
to the dynamics which yields phase space structures to within a given desired accuracy. 

For h — h* sufficiently small and positive, locally the {2n — 1 dimensional) energy surface 
S/i has the structure of 5"^""^ x M in the 2n-dimensional phase space. The energy surface 
S/i is split locally into two components, "reactants" (R) and "products" (P), by a (2n — 2 
dimensional) "dividing surface" (DS(/i)) that is diffeomorphic to S"^""^. The dividing surface 
that we construct has the following properties: 

• The only way that trajectories can evolve from reactants (R) to products (P) (and 
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vice- versa), without leaving the local region £, is through DS(/i). In other words, 
initial conditions on this dividing surface specify all reacting trajectories. 



• The dividing surface is free of local re-crossings; any trajectory which crosses it must 
leave the neighbourhood L before it might possibly cross again. 

• The dividing surface minimizes the (directional) flux. 

The fundamental phase space building block that allows the construction of a dividing 
surface with these properties is a particular Normally Hyperbolic Invariant Manifold (NHIM) 
which, for fixed energy h > h*, will be denoted NHIM(/i). The NHIM(/i) is diffeomorphic to 
S*^""^ and forms the natural dynamical equator of the dividing surface: The dividing surface 
is split by this equator into (2n — 2)-dimensional hemispheres, each diffeomorphic to the open 
(2n — 2)-ball, We will denote these hemispheres by DSf(/i) and DSb(/i) and call them 

the "forward reactive" and "backward reactive" hemispheres, respectively. DSf(/i) is crossed 
by trajectories representing "forward" reactions (from reactants to products), while DSb(/i) 
is crossed by trajectories representing "backward" reactions (from products to reactants). 

The (2n — 3)-dimensional NHIM(/i) is an (unstable) invariant subsystem which, in chem- 



istry terminology, corresponds to the energy surface of the "activated complex" |40l. |59|. 

The NHIM(/i) is of saddle stability type, having (2n — 2)-dimensional stable and unstable 
manifolds IV^lh) and W^{h) that are diffeomorphic to S'^'"'^^ x M. Being of co-dimension 
one {gqI with respect to the energy surface, these invariant manifolds act as separatrices, 
partitioning the energy surface into "reacting" and "non-reacting" parts. 

These phase space structures can be computed via an algorithmic procedure based on 



Poincare-Birkhoff normalization 



38 



39 



52| . This involves developing a new set of coordi- 
nates, the normal form coordinates, {q,p), which are realized through a symplectic coordinate 
transformation from the original, physical coordinates, 

T{q,p) = {q,p), (11) 

which, in a local neighbourhood C of the equilibrium point, "unfolds" the dynamics into a 
"reaction coordinate" and "bath modes". Expressing H in the new coordinates, {q,p), via 

HMq,P) = H{T-\q,p)), (12) 
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gives i^NF in a simplified form. Tfie normalization procedure can also be adapted to yield ex- 
plicit expressions for the coordinate transformations, T{q,p) = {q,p) a.iadT~^{q,p) = {q,p), 
between the normal form (NF) coordinates and the original coordinates [61]. These coordi- 
nate transformations are essential for physical interpretation of the phase space structures 
that we construct in normal form coordinates since they allow us to transform these struc- 
tures back into the original "physical" coordinates. 

The nonresonance condition implies that the normal form procedure yields an explicit 
expression for the normalised Hamiltonian i^^p as a function of n local integrals of motion: 

Hnf = H^F{h, h, ■ ■ ■ , In)- (13) 
The integral, /i, corresponds to a "reaction coordinate" (saddle-type DoF): 

h = qiPi, (14) 
The integrals Ik, for k = 2, ... ,n, correspond to "bath modes" (centre- type DoFs): 

h = l{ql+Pl)- (15) 

The integrals provide a natural definition of the term "mode" that is appropriate in the 
context of reaction, and their existence is a consequence of the (local) integrability in a 
neighborhood of the equilibrium point of saddle-centre-. . .-centre stability type. Moreover, 
the expression of the normal form Hamiltonian in terms of the integrals provides us a way 



to partition the energy between the different modes [62 1. 

The n integrals, the normalized Hamiltonian expressed as a function of the integrals, 
and the transformation between the normal form coordinates and the physical coordinates 
are the key to practically realizing the phase space structures described at the beginning of 
this section. The approximate integrability of Hamilton's equations in the reaction region 
allows a precise and quantitative understanding of all possible trajectories in this region. 
It also provides a natural construction of an energy dependent reaction coordinate whose 
properties are determined solely by the Hamiltonian dynamics, as opposed to the need for 



a priori definitions of possible candidates for reaction coordinates [63 1. 

The n integrals of the motion defined in the neighborhood of the reaction region give rise 
to further phase space structures, and therefore constraints on the motion, in addition to 
those described at the beginning of this section. The common level sets of all the integrals 



are examples of invariant Lagrangian submanifolds GJ, l65|, |66|] , which have the geometrical 
structure of two disjoint n-dimensional toroidal cylinders, denoted MxT"^^, i.e. the cartesian 
product of a line with n — 1 copies of the circle. 

In the next section we show how all of this phase space structure exists for thermostatted 
dynamics of the physical Hamiltonian ([T]) in physical time t. 



IV. MICROCANONICAL PHASE SPACE STRUCTURE: HAMILTONIAN AND 
CORRESPONDING NON-HAMILTONIAN THERMOSTATTED SYSTEMS 

In this section we will show that if the phase space structure described in Section HIT] 
exists for the physical Hamiltonian system ([2]), it also exists in the phase space of the 
non-Hamiltonian isokinetic equations of motion (fTOj) corresponding to the thermostatted 
dynamics of the physical Hamiltonian ([T]) in physical time t. This is accomplished in 3 steps 
by showing: 

1. If the physical Hamiltonian system ([2]) has an equilibrium point at the origin of saddle- 
centre-. . .-centre stability type, then the Hamiltonian system defined by ([5]) corre- 
sponding to the Hamiltonian isokinetic thermostat has an equilibrium point at the 
origin also of saddle-centre-. . .-centre stability type. Morever (and significantly), we 
show that the equilibrium points in these two systems satisfy the same non-resonance 
condition. 

2. The energy of the saddle-centre-. . .-centre type equilibrium point of is negative, 
but it can be brought sufficiently close to zero so that the microcanonical geometrical 
structures described in section HTTl exist on the zero energy surface of ([5]). 

3. The geometrical structures on the zero energy surface of (jS]) map to geometrical struc- 
tures in the phase space of the non-Hamiltonian thermostatted system corresponding 
to (Uni). 

We begin with step 1. We assume that ([2]) has an equilibrium point at (g,p) = {q*,p*) = 
(0, 0). From ([1]), the energy of this equilibrium point is H{0, 0) = $(0). 

The stability of the equilibrium point is determined by the eigenvalues of the derivative 
of the Hamiltonian vector field evaluated at the equilibrium point. This is given by the 
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2n X 2n matrix: 

M... = ^ ^ I , (16) 

y-$gg(0) 0„xr 

where Onxn denotes the nxn matrix of zeros and id„xn denotes the nxn identity matrix. 
We require the equihbrium point to be of saddle-centre-. . .-centre stabihty type. This means 
that the 2n x 2n matrix M^y^ has eigenvalues ±A, ±iuji, i = 2, . . . ,n where A and Ui are real. 

Eigenvalues 7 of M^y, are obtained by solving the characteristic equation det(M3y^ — 
7id2nx2n) = 0. From Theorem 3 of G^l, the block structure of the 2n x 2n matrix M^y^ 
implies that 

det(M,y, - 7id2nx2n) = det(<l>g,(0) + 7'idnxn) = (17) 

so that the 2n eigenvalues 7 are given in terms of a, the eigenvalues of the nxn Hessian 
matrix $gg(0) associated with the potential as follows: 

lk,lk+n = ±V-o-k, k = l,...,n. (18) 

Therefore, if $(g) has a rank-one saddle at g = 0, so that one eigenvalue is strictly nega- 
tive and the rest are strictly positive, then {q,p) = (0,0) is a saddle-centre-. . .-centre type 
equilibrium point for ([2]) as described above. 

We discuss the non-resonance condition in more detail. Suppose ai < and ai > 0. 
i = 2, . . . ,n. T hen the non-resonance condition satisfied by the purely imaginary eigenvalues 
is given by {m,2, . . . , m„) ■ (72, . . . , 7„) 7^ for all integer vectors (m2, . • . , m„) whose entries 
are not all zero (where "■" denotes the scalar product). The non-resonance condition is 
responsible for the existence of the n — 1 (local) integrals of motion /2, . . . , /„. 

Next, we consider the Hamiltonian system ([5]) corresponding to the Hamiltonian isokinetic 
thermostat. It is easy to verify that [q, it) = (0, 0) is an equilibrium point for Q with energy 
/C(0,0) = -|e-*W. 

Proceeding as above, we determine stability of this equilibrium point by computing the 
matrix associated with the linearization of ([5]) at the equilibrium point: 

l^-$,,(0)(ie-*W) 0.x„ ; 

The 2n eigenvalues of Mthcrm, which we denote as 7, can be computed by exactly the same 
type of calculations as above. The resulting eigenvalues are given in terms of the eigenvalues 

11 



of the potential Hessian as follows: 



lk,lk+n = 



iff,* 



1,... 



n. 



(20) 



Therefore, it is clear that if the potential of the physical Hamiltonian, has a rank-one 

saddle at g = 0, so that one eigenvalue is strictly negative and the rest are strictly positive, 
then (g, tt) = (0,0) is a saddle-centre-. . .-centre type equilibrium point for ([5]). 

Moreover, since 7 = ^7 it follows by comparing (fTSll with ( l20l) that if the imaginary 
parts of the eigenvalues associated with the saddle for the physical Hamiltonian satisfy a non- 
resonance condition, then they satisfy a non-resonance condition for the saddle associated 
with the Hamiltonian isokinetic thermostat, i.e (m2, . . . , nin) ■ (72, • ■ ■ , 7n) 7^ implies that 
(ma, ... , m„) ■ ^(72, • • • , 7n) = ("^2, • • • , rrin) ■ (72, . . . , 7„) 7^ 0. 

Now consider step 2. As we showed above, the saddle-centre-. . .-centre type equilibrium 
point (g, tt) = (0,0) of has energy /C(0,0) = —^e'"^^^^ < 0. However, we are only 
interested in the dynamics on the /C = energy surface. The point here is that all of 
the phase space structures described in section UTTl exist for energies "above and sufficiently 
close" to the energy of the saddle-centre-. . .-centre type equilibrium point, and the question 
is whether or not — ie~*^°'' < is close enough to zero so that the phase space structures 
described in Section UTTl exist on the /C = energy surface. This can easily be arranged by 
making $(0) larger by adding an appropriate constant to $(g), or by changing the value of 
the temperature T. 

The final step 3 is to show that the phase space structure of ([5]) on /C = exists for 
the isokinetic equations of motion ffTOj) corresponding to the thermostatted dynamics of the 
physical Hamiltonian ([T]) in physical time t. 

Step 3 follows from general results that show that invariant manifolds and their sta- 
bility properties are preserved under differentiable, invertible (with differentiable inverse) 
coordinate transformations (i.e. they are preserved under differentiable conjugacies). We 
emphasize that these results are "well-known" and appear in a variety of places throughout 



These general results allow us to make the following conclusions: 

• Under the map ([7]) the 2n — 1 dimensional invariant energy surface of the effective 
Hamiltonian system ([5]) maps to a 2n — 1 dimensional invariant manifold for the non- 



the literature, e.g., j68|. 
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Hamiltonian isokinetic equations of motion for the thermostatted physical Hamiltonian 
([1]) in physical time t defined by ffTOj) . 

• Under the map ([7]) the 2n — 3 dimensional NHIM, its 2n — 2 dimensional stable and 
unstable manifolds, the n-dimensional invariant Lagrangian submanifolds, and the 
2n — 2 dimensional dividing surface map to a 2n — 3 dimensional NHIM, its 2n — 2 
dimensional stable and unstable manifolds, n-dimensional invariant submanifolds, and 
a 2ri — 2 dimensional dividing surface in the 2n — 1 dimensional invariant manifold in 
the 2n dimensional phase space of the non-Hamiltonian isokinetic equations of motion 
for the thermostatted physical Hamiltonian ([1]) in physical time t defined by ffTOl) . 



V. SUMMARY AND OUTLOOK 



In this paper we have examined the relation between phase space structures in Hamil- 
tonian and non-Hamiltonian thermostats. In particular, we have established the existence 
of a mapping between invariant phase space structures in the phase space of the extended 
Hamiltonian for the isokinetic thermostat and corresponding structures in the phase space 
of the non-Hamiltonian Gaussian isokinetic thermostat. 

Our results establish a conceptual link between the question of thermostat ergodicity and 
the issue of statisticality in unimolecular isomerization reactions. The existence of normally 
hyperbolic invariant manifolds in both the physical and extended Hamiltonian phase spaces 
means that recently developed methods for the analysis of isomerization dynamics can be 
applied to the thermostat problem. Numerical studies based on the ideas presented here are 
currently in progress. 



Finally, we note that the approach presented here s 
tonian thermostats, such as the Nose-Poincare system 



lould be applicable to other Hamil- 
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